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1  Introduction 


Near-field  (NF)  antenna  measurement  techniques  have  been  an  indispensable  tool  for  deter¬ 
mining  antenna  far-held  (FF)  patterns  over  many  decades  [1]  [2] .  NF  data  are  almost  always 
collected  on  a  planar,  cylindrical,  or  spherical  surface  mainly  to  facilitate  FF  computation. 
The  cylindrical  and  spherical  scanning  techniques  are  often  used  to  determine  FF  patterns 
over  a  wide  range  of  directions,  while  the  planar  technique  is  reserved  for  antennas  with 
narrow  FF  patterns.  NF  techniques  have  also  been  explored  to  determine  the  bistatic  RCS 
of  a  target.  Cote  and  Wing  [3]  demonstrated  that  the  bistatic  RCS  can  be  computed  from 
spherical  NF  scan  data.  More  recently,  Hansen  et  al.  [4]  extended  for  target  scattering  the 
standard  cylindrical  NF  theory  that  was  originally  developed  for  antennas,  while  Marr  et 
al.  [5]  experimentally  demonstrated  the  feasibility  of  determining  the  bistatic  RCS  from 
cylindrical  NF  data.  Cowan  and  Ryan  [6]  and  Farr  et  al.  [7]  applied  the  planar  NF  scanning 
technique  to  bistatic  scattering  from  targets,  while  Zahn  and  Sarabandi  [8]  applied  it  to 
bistatic  scattering  from  random  rough  surfaces. 

One  significant  difference  between  antenna  and  target  NF  measurements  is  that  antennas 
are  excited  by  internal  sources,  while  targets  are  excited  by  a  plane  wave  that  is  generated  by 
a  source  external  to  them.  In  target  NF  measurements,  the  NF  probe,  therefore,  cannot  be 
placed  in  the  region  between  the  target  and  the  transmit  dish  that  generates  the  plane  wave, 
as  its  placement  would  distort  the  incident  plane  wave.  Hence,  for  bistatic  RCS  applications, 
it  is  not  possible  to  collect  scattered  NF  samples  over  the  entire  47T  steradian  in  the  spherical 
NF  scanning,  or  over  [0,  2n)  in  the  cylindrical  scanning.  Since  only  truncated  NF  samples 
are  available,  computed  FF  solutions  suffer  from  the  angular  truncation  error;  if  NF  samples 
are  collected  over  [0i,  02]  on  a  scan  surface  of  radius  a,  then  FF  values  can  be  computed 
accurately  only  over  [03,04]  where  0 3  >  0i  and  04  <  02.  Hansen  et  al.  [4]  showed  that  the 
extent  of  the  truncation  error,  as  determined  by  A0i  =  03  —  0i  and  A02  =  02  —  04,  depends 
on  the  scan  radius. 

Antenna  FF  patterns  computed  from  truncated  NF  data  are  subject  to  similar  trunca¬ 
tion  errors.  With  antennas,  it  is  possible,  if  desired,  to  completely  eliminate  the  angular 
truncation  error.  With  the  spherical  scanning,  radiation  patterns  in  all  directions  can  be 
determined,  if  NF  samples  are  collected  at  an  appropriate  sampling  rate  over  the  entire 
surface  of  a  sphere.  Similarly,  with  the  cylindrical  scanning,  patterns  can  be  determined  in 
all  azimuthal  (horizontal)  directions  over  a  finite  vertical  extent,  if  NF  data  are  collected 
over  the  entire  vertical  surface  of  a  finite- height,  right-circular  cylinder.  Since  FF  angular- 
truncation  errors  in  spherical  and  cylindrical  scans  can  be  eliminated  (albeit  at  the  expense 
of  substantially  increased  data  collection  time),  it  appears  that  much  of  the  existing  body 
of  research  on  truncation-error  reduction  has  been  concerned  with  reducing  FF  errors  that 
result  when  the  linear  extent  of  a  scan  surface  is  truncated,  as  occurs  in  the  planar  scanning 
or  the  height  truncation  in  the  cylindrical  scanning  [4]  [10]  [11]  [12]  [13]. 

I11  this  paper,  we  consider  angular-truncation  error  reduction  for  the  2D  cylindrical/spherical 
NF  scanning.  After  examining  the  limitations  of  the  traditional  multipole-based  expansion 
of  angular-truncated  scan  data,  we  consider  an  alternative  expansion  based  on  Slepian  func- 
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tions  [14]  and  develop  a  NF-to-FF  transformation  for  the  resulting  expansion  coefficients. 
We  apply  this  new  NF-to-FF  transformation  to  2D  bistatic  scattering  from  a  homogeneous 
dielectric  circular  cylinder  and  establish  that,  compared  with  the  traditional  multipole-based 
NF-to-FF  transformation,  it  significantly  reduces  the  angular-truncation  error. 
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2  Truncation  Error  in  the  Traditional  NF-to-FF  Trans¬ 
formation 

In  the  traditional  2D  cylindrical/spherical  NF  scanning,  the  scan  data,  E(a,(ft),  collected 
over  a  circle  of  radius,  a,  is  expanded  in  terms  of  circular  multipoles,  Hm(ka)e~im^ ,  where 
Hm(ka )  is  the  mth-order  Hankel  function  of  the  first  kind  and  k  is  the  wavenumber: 

M  M 

E(aA)=YJ  bmHm(ka)e~im^=J2  cm(a)  e"im*  (1) 

m=—M  m=—M 

Here,  M  =  int(fca0)  +  N0,  where  int  is  the  integer  function,  a0  is  the  radius  of  the  minimum 

circle  that  completely  encloses  the  antenna/scatterer,  and  N0  is  a  small  integer.  If  the  scan 
data  is  available  over  the  entire  interval,  [0,  27t),  then  the  expansion  coefficients,  bm,  can  be 
determined 

r2  7r 

bm  =  O  TT  a.  \  /  E ^  -  M<m<M ,  (2) 

ZTrHm{ka)  Jo 

and  the  radiated/scattered  held,  E(p,(ft),  can  be  computed  anywhere  for  any  p  >  a0  and 
(ft.  If  the  scan  data  is  collected  over  a  limited  angular  region,  [(ft0,  27t  —  (ft0],  which  is  often 
the  case,  it  is  customary  to  assume  the  NF  values  outside  the  angular  region  to  be  zero  and 
obtain  the  expansion  coefficient,  brn.  from  the  truncated  NF  scan  data 

bm  =  '  r^°  E(a,  (ft)eim*d(ft,  —  M  <m  <  M.  (3) 

2nHm(ka)  J<$,0 

Since  bm  ^  bm:  the  FF,  lim  E(p,<f)),  computed  from  bm  is  not  expected  to  agree  with  the 
FF,  lim  E(p,(f>),  computed  from  bm  over  the  entire  angular  region,  [0,  27t). 

Using  (1)  and  (3),  we  may  identify  the  angular  region  where  lim  E(p,  (ft)  and  lim  E(p,  (ft) 
are  expected  to  disagree: 

M 

lim  E(p,  <ft)  =  liin  bmHm(kp)e~im(t> 

P^°°  P^°°  m=—M 

1  I  2  /*27T — (f>0 

=  ekt  (4) 

M 

Here,  Tm(cl,  (ft— (ft')  =  ^  _7r/2) / Hm(ka),  which  is  referred  to  as  the  taper  function  in 

m=—M 

[4],  determines  the  angular  extent  of  NF  scan  data  E(a,  (ft1)  that  is  required  to  compute  the 
FF,  lim  E(p,  (ft).  Figure  1  plots  normalized  | Tm{cl,  (ft  —  (ft')\  as  a  function  of  (ft— (ft1  for  various 

values  of  scan  radius,  a  and  M  =  90.  For  a  =  1000A,  \Tm(o>,  (ft  —  <ft')\  is  sharply  peaked  at 
(ft  =  eft'.  However,  as  a  gets  smaller,  |Tjvf(a,  (ft— (ft')\  develops  a  broad  plateau  around  (ft  =  (ft'. 
Thus,  as  a  gets  smaller,  more  extended  NF  scan  data  around  (ft  =  <ft'  is  needed  to  compute 


3 

Approved  for  public  release;  distribution  unlimited. 


-I5D  -IDD  -5D  D  5D  1DD  15D 


<j>  -  <t>  (deg) 


Figure  1:  Normalized  taper  function  |Tm(o,  0  —  <p')\  for  various  values  of  scan  radius,  a  and  M=90. 

lim  E(p,  0).  This  behavior  of  \Tm(cl,  0—0') |  for  small  a  is  responsible  for  the  truncation  error 

in  the  FF.  Even  though  the  width  of  the  plateau  near  0  =  0'  decreases  as  a  increases,  thus 
reducing  the  FF  truncation  error,  it  generally  takes  longer  to  collect  NF  data  at  a  larger  scan 
radius  (The  number  of  required  NF  samples  is  independent  of  a.  However,  the  NF  probe 
must  traverse  longer  distances  at  large  a)  . 
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3  Alternative  Expansion  of  NF  Scan  Data 

We  may  observe  that  the  fundamental  reason  for  the  FF  truncation  error  is  that  the  multipole 
field,  Hm(ka)e~im^ ,  that  is  used  to  expand  the  NF  scan  data,  fails  to  be  orthogonal  over 
[4>o,  2i T  —  (f>0\,  leading  to  bm(a)  ^  bm(a).  Therefore,  we  may  seek  an  alternative  expansion  of 
the  scan  data  in  terms  of  a  basis  function  set  that  is  orthogonal  over  [</>„,  27 T—(p0].  It  is  well 
known  that  the  radiated/scattered  field  is  spatially  band-limited  [9].  Stated  in  the  context 
of  (1),  it  is  index-limited,  i.e.,  cm(a)  =  0  for  \m\  »  M.  When  the  expansion  coefficients, 
cm(a),  are  Fourier-transformed  to  angular  domain,  using  (1),  to  obtain  the  radiated/scattered 
field,  the  resulting  field  is  not  angle-limited,  i.e.,  \E(a,4>)\  >0,  0  <  (ft  <  2i r,  thus  making 
it  difficult  to  find  a  basis  function  set  that  is  orthogonal  over  [<p0,  2tt  —  <p0\  and  satisfies  the 
Helmholtz  equation. 

It  is  possible  to  find  a  basis  function  set  that  results  in  the  maximum  concentration  of 
’’energy”  in  a  given  truncated  angular  domain.  Following  [14],  we  define 

1  r  27r  <f>o 

S (a,  <f>0)  =  —  E*  (a,  0)  E(a,  (p)  d<p,  (5) 

Zll  J  (fi0 

and  seek  a  basis  function  set  for  E(a,(p)  that  maximizes  S(a,cp0 )  for  given  (pa.  Using  (1), 

S(a,  (t>0)  =  £  4(a)£iW(0o)  <w(a)  =  C+ -T -C  >  0,  (6) 

m  m! 

where  C  is  the  column  vector  containing  the  2 M  +  1  values  of  cm(a);  C+ ,  its  hermitian 
conjugate;  and  K ,  the  (2 M  +  1)  x  (2 M  +  1)  real,  symmetric,  Toeplitz  matrix  with 

Km,m,  =  —  /  e<m~m  *dxj>  =  -  ^  sine  (m  -  m')<j>0,  (7) 

ATT  J (j)0  7T 

— M  <  m,ml  <  M,  and  sinc(x)  =  sin(a:)/a;.  Since  S(a,  (p0)  is  always  positive,  K  is  a  positive- 
definite  matrix  [14].  Let  vn  and  An  be  the  nth  eigenfunction  and  eigenvalue  of  K  so  that 
K-vn  =  Xnvn,  n  —  1,  ...  ,2M+1.  Or 


K  •  V  =  V  ■  A, 


(8) 


where  V  =  [ui  v2  v3  ...  v2m-i  v2m  ^2M+i]  and  A  is  the  diagonal  matrix  with  Ann  =  Xn, 
n  —  1,  ...  ,  2M+1.  Since  K  is  positive-definite,  V  is  a  real  matrix  with  V  —  V  where  V 


is  the  transpose  of  V,  and  Xn  is  real  with  1  >  Ai  >  A2  >  ...  >  X2m  >  X2m+i  >  0  [14]. 

We  may  use  vn  to  construct  an  expansion  function  set  that  is  orthogonal  over  [(p0,  2i T—(p0\. 
We  define 

M 

Y,  vn(m)  ,  n  =  1,  ....,  2M+1  (9) 

m=—M 
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(A):  eigenvalue  distribution 

50 1 - ! - 


(B):  orthogonality  of  sn  (§) 


0 


Figure  2:  (A)  distribution  of  eigenvalues  (B)  orthogonality  of  sn(0) 


where  vn(m)  is  the  mth  element  of  vn .  sn(4>)  is  commonly  called  the  prolate  spheroidal  wave 


function, 

or  Slepian  function  [14].  It  satisfies  [14] 

r27 T—(f>o 

1  ^n(^)  ^n'  ($)(%(/)  ^n^n,n' 

J  (j)0 

(10) 

and 

p  2*7T 

/  S*n(4>)  Sn/(0)d0  =  Sny 

Jo 

(11) 

According  to  (10),  the  eigenvalue,  An,  physically  corresponds  to  the  concentration  of  sn(4>) 
over  [0o,27t  -  (p0\  [14]. 

Plotted  in  Figure  2. A  is  the  distribution  of  eigenvalues  corresponding  to  M  =  90  and 
4>0  =  30°,  where  we  have  introduced  a  floor  of  -156  dB  to  avoid  having  to  compute  the 
logarithm  of  0.  If  the  sn(0)  are  chosen  to  maximize  fzl  E(a,  (j))E*(a,  (j>)d(j),  then  there  are 
approximately  2M(j)0/n  eigenvalues  that  are  close  to  1  [14].  Since  we  have  instead  chosen  the 
sn(<j))  to  maximize  J^~^°  E(a,  (j))E*(a,  4>)d(ft,  we  expect  approximately  2M(1  —  4>0/ti )  =  150 
eigenvalues  close  to  1,  as  shown  in  the  figure.  Figure  2.B  shows  the  orthogonality  of  the 

sn(4>). 
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Since  the  sn(cft)  are  orthogonal  over  [cft0,  2i T  —  (ft0],  we  may  use  them  to  expand  the  NF 
scan  data,  E(a,  (ft): 

p 

E{a,  </>)  =  5Z  d„(a)  s«(^)  —  ^(a)  '  s(</>)  (12) 

n=l 

where  P  <  2 M  +  1;  d(a)  and  s(<ft)  are  the  column  vectors  of  length  P  containing  the  dn(a) 
and  sn((ft)  values,  respectively.  Using  (10),  the  expansion  coefficients,  d(a),  can  be  obtained 
from  the  truncated  NF  scan  data: 

r2-K-<f>0 

dn(a)  =  —  E(a,cft)s*n((ft)d(ft ,  n  =  1,  ...  ,P.  (13) 

An  J  <fio 

Since  An  — >  O"1"  as  n  — >  2M+1,  P  is  chosen  so  that  A p  >  e,  where  e  is  a  small  positive  number. 
We  note  that  dn(a)  can  be  expressed  in  terms  of  the  multipole-expansion  coefficients,  bm. 
Substituting  (1)  and  (9)  into  (13),  and  using  (8), 

d(a)  =  V2 7T  Vp  ■  T(ka)  •  b  (14) 

where  d(a)  and  b  are  the  column  vectors  containing  the  P  values  of  dn(a)  and  2 M  +  1  values 
of  bm,  respectively;  Up  is  the  (2 M  +  1)  x  P  matrix  obtained  by  taking  the  hrst  P  columns 
of  U;  and  H(ka )  is  the  (2 M  +  1)  x  (2 M  +  1)  diagonal  matrix  whose  diagonal  elements  are 
Hm(ka)  with  —  M  <  m  <  M.  With  P  =  2 M  +  1  in  (14),  it  is  evident  that  dn(a)  remains 
finite  for  all  n,  in  spite  of  Xn  — »  0+  as  n  — >  2M+1. 

Substitution  of  (14)  into  (12)  yields 

E(a,  (ft)  =  J2  bmHm(ka )  e~im'+  (15) 

m  m' 

P 

where  [Qp\m,m'  =  X]  vn(m )  vn(m!).  We  note  that  if  P  =  2M  +  1,  [Qp]mim/  =  and  thus 

n=l 

(12)  satishes  the  2D  Helmholtz  equation  rigorously.  If  P  <  2 M  +  1,  then  [QP\m,m'  ^ 
but  has  a  strong  diagonal  dominance  as  shown  in  Figure  3  for  P  =  156  and  166  for  the  sn(cft) 
considered  above.  Therefore,  when  P  <  2 M  +  1,  (12)  satishes  the  2D  Helmholtz  equation 
only  approximately. 
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4  Transformation  of  the  Expansion  Coefficients,  dn(a) 

In  order  to  obtain  FF  from  the  NF  representation,  (12),  we  define  the  operator.  T (p,  a),  that 
propagates  the  scan  data,  E(a,  (ft),  to  the  field  at  (p,  (ft): 

p 

E(p,  (ft)  =  T(p,  a)  ■  E(a ,  (ft)  =  J2  [T(/h  a)  dn(a )]  sn((ft) 

n=  1 

=  51  ddp)  sn{(ft)  =  dT (p)  ■  s((ft)  (16) 

n= 1 


Using  V  P  -V  p  =  Ip,  where  Ip  is  the  identity  matrix  of  order  P,  we  establish  that  Up  is  the 

=  T 

pseudo-inverse  of  V P  [15]  and  obtain  the  least-square  solution  of  (14), 


bLS  = -L=H  {ka)  ■  V P  ■  d(a) . 
v  27T 

Substitution  of  the  above  equation  into  (14)  yields 


d(p)  =  T(p,  a)  ■  d(a)  = 


V P  ■  H(p,a)  ■  V P 


d(a), 


(17) 


(18) 


where  T(p,a)  is  the  matrix  representation  of  T(p,  a);  and  H(p,a)  is  the  (2M  + l)x(2M  + 
1)  diagonal  matrix  whose  diagonal  elements  are  Hm{kp) / Hm(ka),  —M  <  m  <  M.  The 
FF  expansion  coefficients,  lim  dn(p),  are  then  obtained  by  substituting  the  large-argument 

expression  for  Hin(kp),  or 


lim  dip) 

p—>  OO 


Up  ■  Goo  (a)  ■  Up 


■  d{a), 


(19) 


where  G^a)  is  the  (2Md)x(2Md)  diagonal  matrix  whose  diagonal  elements  are  e  / Hm{ka), 
—M <m<  M.  We  note  that  as  the  scan  radius,  a,  approaches  infinity,  lim  Hip,  a)  =  Pm+i, 

p,a—>  oo 

where  Hm+i  is  the  identity  matrix  of  order  2 M  +  1.  Since  VP  -  Up  =  Ip,  we  have 
lim  d(a)  =  limd(p)  and  consequently  lim  E(a,(ft)  =  lim  E(p,(ft),  satisfying  the  basic  physi¬ 
cal  requirement  that  the  scan  data  collected  at  a  =  00  correspond  to  the  FF. 

In  the  traditional  2D  NF-to-FF  transformation  based  on  (3),  both  the  the  expansion 
coefficients,  bm,  and  the  FF  can  be  computed  using  0(Mlog2M)  floating-point  operations 
(FPOs)  by  taking  advantage  of  the  Fast  Fourier  Transformation  (FFT).  The  new  NF-to- 
FF  transformation  algorithm,  however,  requires  the  eigen-decomposition  of  K,  which  re¬ 
quires  0(M3)  FPOs.  Thus,  from  the  purely  computational  efficiency  point  of  view,  the  new 
NF-to-FF  transformation  algorithm  is  less  attractive  than  the  traditional  multipolc-based 
algorithm. 
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5  Numerical  Examples 

We  apply  the  new  NF-to-FF  transformation  to  the  2D  bistatic  scattering  scenario  shown  in 
Fig.  4.  A  plane  wave  of  wavelength,  A,  propagates  to  the  left  and  is  scattered  by  a  circular 
cylinder  of  radius  aD  =  10A  and  dielectric  constant  1.6.  The  incident  electric  held  is  polarized 
perpendicular  to  the  cylinder  axis.  The  exact  bistatic  near  and  far  fields  can  be  computed 
using  (1)  with  the  expansion  coefficients,  bm,  given  in  [16]  with  M  =  90.  To  compute  far 
fields  from  near  fields,  we  sample  the  scattered  NF,  P(a,0),  using  an  ideal  probe  on  the 
concentric  cylinder  of  radius  a  =  25A  from  0  =  30°  to  0  =  330°  with  increments  of  A0  =  O.5°, 
which  is  about  one  fourth  of  the  Nyquist  sampling  interval,  A 0  =  7t/(M  +  1)  [1] [2] .  We 
compute  two  sets  of  FF  values  from  the  NF  samples  using  (3)  and  (19),  respectively.  (3)  is 
approximated  by 


A0 


j 


2n Hm(ka) 


E(a,  0j)e<m^ ,  -  M  <  m  <  M, 


where  J  is  the  number  of  NF  samples  and  the  summation  over  j  can  be  evaluated  efficiently 
using  FFT.  To  compute  the  FF  using  (19),  we  first  construct  the  181x181  K  matrix  taking 
advantage  of  its  symmetric  Toeplitz  structure.  Since  the  integral  in  (13)  needs  to  be  evaluated 
numerically,  the  matrix  element,  Km  mi ,  is  also  evaluated  numerically  using  the  Trapezoidal 
rule  rather  than  using  the  analytic  expression  given  in  (7).  We  then  compute  the  eigenvalues, 
An,  and  eigenvectors  vn,n  =  1,  ...  ,  181.  The  resulting  Xn  distribution,  as  shown  in  Fig.  2. A, 
contains  P=165  non-negligible  eigenvalues  with  Xn  >  10~14.  Thus,  we  compute  the  Slepian 
functions,  sn(0),  and  the  coefficients,  dn(a ),  for  n  =  1,  ...,  165  to  expand  the  NF  scan  data. 
The  FF  expansion  coefficients,  lim  dn(p),  and  FF  values,  lim£'(p,  0),  are  obtained  using 

(19)  and  (16),  respectively. 

The  three  sets  of  normalized  FF  solutions  are  plotted  in  Fig.  5  as  a  function  of  bistatic 
angle.  As  expected  from  the  behavior  of  the  taper  function,  T(a,  0  —  00,  the  FF  solution 
computed  using  the  traditional  NF-to-FF  transformation  deviates  from  the  reference  solu¬ 
tion  for  30°  <  0  <  62°.  In  contrast,  the  FF  solution  computed  using  the  new  NF-to-FF 
transformation  diverges  from  the  reference  solution  only  for  30°  <  0  <  33°,  demonstrating 
that  the  new  NF-to-FF  transformation  significantly  reduces  the  error  in  the  computed  FF. 
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30° 


a 


330° 


Figure  4:  2D  NF  scanning  for  bistatic  scattering  applications. 


It  is  well  known  that  the  performance  of  many  numerical  techniques  that  take  advantage 
of  the  eigen-structure  of  a  matrix  depends  on  signal-to- noise  ratios  (SNR).  One  significant 


advantage  of  NF  measurements  is  that  it  is  possible  to  achieve  a  high  SNR  [2] ;  SNRs  of  40 
to  60  dB  are  routinely  achieved  [17].  In  order  to  investigate  the  dependence  of  performance 
on  SNR,  we  introduce  noise  to  the  exact  NF  scan  data,  E0(a,(pj ): 


E(a,  (j>j) 


E0(a,  c t>j )  +  N. 


Tji  j  1 )  ...,</ 


where  Nj  is  modeled  as  a  complex,  zero- mean,  Gaussian  random  process.  Compared  in  Fig. 
6  are  the  two  normalized  FF  solutions  computed,  respectively,  using  the  new  and  traditional 
NF-to-FF  algorithms  for  the  scattering  problem  of  Fig.  4  with  SNR=43  dB.  The  figure  shows 
that  the  addition  of  noise  degrades  the  performance  of  the  new  algorithm  more  significantly 
than  that  of  the  traditional  algorithm;  the  solution  generated  with  the  new  algorithm  now 
agrees  with  the  reference  solution  over  </>> 45°,  while  the  noise-free  solution  shown  in  Figure 
5  agrees  with  the  reference  solution  0>  33°.  Even  with  its  higher  sensitivity  to  noise,  the  new 
NF-to-FF  algorithm  still  produces  the  FF  solution  that  agrees  with  the  reference  solution 
over  a  wider  range  of  bistatic  angle. 

If  we  let  E(p)  and  E(a)  be  the  column  vectors  containing  the  J  values  of  E(p,  0)  and 
E(a,4>),  respectively,  the  traditional  algorithm,  (3),  may  be  expressed  in  matrix  form  as 


(20) 


and  the  new  algorithm,  (16),  as 


(21) 
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Figure  5:  FF  solutions  computed  from  NF  scan  data.  Top  figure:  Normalized  FF  as  a  function  of  bistatic 
angle.  Bottom  figure:  Same  as  the  top  figure  except  that  the  FF  solutions  are  plotted  from  30  to  80°. 
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Figure  6:  FF  solutions  generated  from  NF  scan  data  with  SNR=43dB. 
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Figure  7:  Comparison  of  the  transformation  matrices:  (A)  101og10  rtr|,  (B)  101og10  \  fnew\  with  P=156 
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Here,  F  denotes  the  Fourier  transformation  matrix  of  (9)  and  QP  =  V pV  P,  and  U P  = 
_  — —  —  i _ t  =-i  =  -i 

V pAP  V Pl  where  AP  is  the  diagonal  matrix  of  order  P  with  [AP  ]n^n  =  1/An.  One  may  use 
(21)  to  compute  the  E(p)  without  explicitly  constructing  the  sn(4>).  Figure  7  compares  the 
two  ’’transformation  matrices”  for  the  NF  scan  scenario  considered  above  with  p  =  2.5x10' A. 
The  Ttr  is  Toeplitz,  as  required  by  (4).  The  rnew  is  plotted  for  P  — 156.  For  60°  <  0  <  300°, 
the  Ttr  and  rnew  behave  quite  similarly.  Outside  this  angular  region,  where  the  angular 
truncation  error  is  of  concern,  the  Tnew  attempts  to  ’’collect  maximum  information”  from 
all  available  E(a,4>')  to  construct  the  E(p,cp).  It  is  to  be  noted  that  the  Tnew  shows  a  high 

sensitivity  to  P  for  0<  60°  and  0>3OO°  due  to  the  presence  of  AP  through  U p  in  (21).  This 
explains  the  higher  noise  sensitivity  of  the  FF  solution  obtained  using  the  new  algorithm. 
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6  Conclusion 


We  have  introduced  a  new  NF-to-FF  transformation  algorithm  for  the  2D  cylindrical/spherical 
NF  scanning.  Compared  with  the  traditional  multipole-expansion-based  algorithm,  this  new 
algorithm  significantly  reduces  the  error  in  the  FF  solution  that  results  when  NF  data  is 
truncated,  albeit  at  a  higher  computational  cost.  Even  though  we  were  primarily  motivated 
by  bistatic  RCS  applications,  the  new  algorithm,  when  fully  extended  to  3D,  can  be  used  to 
reduce  the  required  scan  area  in  cylindrical  and  spherical  antenna  NF  measurements,  thereby 
making  these  measurements  more  cost-effective. 
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Acronyms 


NF  Near  Field 

FF  Far  Field 

RCS  Radar  Cross  Section 
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